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Abstract 

In this paper, we formulate a generalization of the simple Bekenstein- 
Sandvik-Barrow-Magueijo (BSBM) theory of varying alpha by allowing 
the coupling constant, oj, for the corresponding scalar field to depend on 
ip. We focus on the situation where cj is exponential in ?/) and find the late- 
time behaviours that occur in matter-dominated and dark-energy domi- 
nated cosmologies. We also consider the situation when the background 
expansion scale factor of the universe evolves in proportion to an arbi- 
trary power of the cosmic time. We find the conditions under which the 
fine structure 'constant' increases with time, as in the BSBM theory, and 
establish a cosmic no-hair behaviour for accelerating universes. We also 
find the conditions under which the fine structure 'constant' can decrease 
with time and compare the whole family of models with astronomical data 
from quasar absorption spectra. 

PACs nos: 98.80.Es, 98.80.Bp, 98.80.Cq 

I. INTRODUCTION 

The electromagnetic fine-structure constant, a — e'^/hc, has traditionally 
been thought of as a fundamental constant of Nature. However, high-redshift 
observations of quasar absorption spectra have continued to suggest [TJ dl 131 HI IS] 
that this 'constant' may exhibit very slow temporal and spatial variation. One 
theory that generalises Maxwell's theory of electromagnetism and the general 
theory of relativity in order to accommodate and test these possible variations 
in a is the BSBM theory of Bekenstein and Sandvik, Barrow and Magueijo 
[51 El [S]. In this theory, space-time variations in a are carried by a scalar 
field which couples to electrically charged matter; this is an example of what 
has since become known as a 'chameleonic' scalar field. These variations, in 
turn, contribute to the space-time geometry and cosmological dynamics. Such a 
theoretical development is important because almost all published observational 
bounds on variations in a derive from simply allowing a to become a variable 
in the conventional equations of physics in which it is strictly a constant. The 



BSBM theory provides a self-consistent varying-a theory in the way that Brans- 
Dicke theory [5] describes a varying Newtonian gravitation 'constant', G, via a 
covariantly conserved Brans-Dicke scalar field (j)BD ^ G^^. 

In the BSBM theory, the quantities c and h are taken to be constants, and 
variations in a are ascribed to changes in e, the electron charge. We write 
e = eoe'^, where tjj = ipix'^) is a dimensionless scalar field, and eo is the present 
value of e. The BSBM Lagrangian is [7] 

>C = ^/^{Lg + Cmat + l^i, + l^erne-^'^'^) (1) 

where Zg = = -^d^tpd^ip and £em = "zUff^"- Here, i? is the Ricci 

curvature scalar, a; is a coupling constant, and we have defined an auxiliary 
gauge potential = with corresponding field tensor /^^ — e^'Ff^i, = 

d^an — d^a^, so that the covariant derivative takes its usual form — 9^ + 
ieoa^. The quantity Cem is usually parametrized by the ratio ^ — Cem/ P, where 
p is the total baryon energy density. This ratio, C, describes the fraction of non- 
relativistic matter in the universe that contributes to £cm- The cosmological 
value of ^ (denoted Cm) has to take into account non-baryonic matter, and thus 
depends strongly on the nature of dark matter, the nature and constituents of 
which are currently uncertain. In general, therefore, Cm is determined by the 
relative role of magnetic, B, and electric, E, fields in the dominant cold dark 
matter content of the universe, and Cm = {E^ — B'^)/{E^ + B^) can take any 
value between —1 and 1 [71,|H]; also, during the radiation era, (Cm) = 0. In 
the situations where Cm < 0, a slow logarithmic increase in a occurs during 
the matter-dominated era but a then tends to a constant when dark energy 
causes the expansion of the universe to accelerate [8] . Therefore, if a is found to 
decrease with time during the matter-dominated era, the simple BSBM model 
with Cm < would be inconsistent with this variation. 

Observational searches for varying a at cosmological redshifts have used 
high-precision studies of quasar spectra coupled with many-body calculations of 
the effects of small changes to a on relativistic corrections to atomic transition 
frequencies [10]. The most recent analysis by Webb et al |5] of a new set of 
observations of quasar spectra from the VLT gives the following best fits for 
Aa/a = {q:{z) — q;o)/q;o, where ao denotes the present value of a and z is the 
redshift of the absorption lines: 

Aa/a = (-0.06 ± 0.16) x lO^^ for z < 1.8, 
Aa/a = (0.61 ± 0.20) x 10"^ for z > 1.8. 

For comparison, an earlier analysis of data from the Keck Telescope gave 

DP- ia> 01: 

Aa/a ^ (-0.54 ± 0.12) x 10"^ for z < 1.8, 
Aa/a = (-0.74 ± 0.17) x 10"^ for z > 1.8. 
All the errors quoted are la errors. Note that, for redshifts z > 1.8, the 
sense in which a{z) varies is different for these two sets of observations. The 
two samples were obtained from telescopes in different hemispheres of the Earth, 
which has led to a proposal in [S] that there may be a strong spatial (dipolar) 
dependence in the variation of a. However, the sense of the overall temporal 
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variation is unclear, and it is conceivable that a could increase, decrease or re- 
main constant with time during the matter era. Further evidence that a may 
increase during part of the matter-dominated era can be found from observa- 
tions based on single ion differential measurements, which place a constraint of 
Aa/a = (0.566±0.267) x IQ-^ at z = 1.84 [TT], [12] that the authors interpret as 
an upper limit because of unknown systematics. For a summary of these results, 
and others, the reader is referred to the recent theoretical reviews by Martins 
[13] and Uzan [M]. However, note that the observational results announced by 
Srianand et al [15], which appeared at first to be consistent with no significant 
redshift variation in a, have been shown to have methodological problems, and 
a subsequent re-analysis of the same data produced results consistent with the 
variations found in [TJ [H |31 13] . For a discussion of some of the observational 
problems, see [16]. 

In this paper, we will show how the simple BSBM theory can be gener- 
alised in a natural manner by allowing the coupling constant lo in to take 
on a V'-dependence, so that uj = We will show that a simple choice 

of the coupling function, uj{ijj), permits Aa/a to increase or decrease during 
the matter-dominated era according to the value of a single parameter which 
vanishes when the new theory reduces to the BSBM case (with constant uj). 
However, there is still a general pattern to the late-time a{t) behaviour in an 
expanding universe, and solutions for a{t) generally tend to a particular asymp- 
totic form which is independent of the initial conditions. We shall focus our 
analysis on the late-time behaviour of the matter-dominated and A-dominated 
eras of a zero-curvature Friedmann universe, but also provide solutions for the 
more general case in which the scale factor of the universe evolves as a{t) = t", 
where n is a constant, which mimics negative curvature dominated expansion 
when n = 1. We find that, as in the original BSBM theory, this extended 
theory still predicts that any variation in a is rapidly suppressed during the A- 
dominated era, thus allowing complex atoms to persist into the far future (see 
[17]). Our generalised theory also gives more scope for spatial variations in a 
to be described because gradients in a can now be driven by gradients in Lu^tp). 

In section II we set up the basic equations of the generalized BSBM theory 
and explore some of their key properties. In section HI we discuss the choice of 
^(-0) and in section IV we study the cosmological solutions in a dust-dominated 
flat Friedmann universe. In section V we compare this with the behaviour in 
a wider range of Friedmann background metrics. In section VI we investigate 
what happens to the variation of a when the universe undergoes accelerated 
expansion, dominated by dark energy in the form of a positive cosmological 
constant. In section VII we summarise these results and compare them with the 
cosmological data. Finally, section VIII contains some concluding discussion. 

II. THE UNDERLYING MODEL AND BACKGROUND 

EQUATIONS 

Consider the modified BSBM Lagrangian 

C = ^/^{Cg + £mat + C^, + Ccme'^"^) (2) 
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where Zg = = -'^d^.tpd'^ip and £om = -Iff.uf'"', where the 

couphng uj is now taken to be a function of the scalar field "0. In order to ensure 
that the energy density of the scalar field, p = ^^^^^■0^, is non-negative, we will 
impose the no-ghost condition uj{ip) ^ 0. As in the BSBM theory, the fine 
structure 'constant' is determined by 

a = Qfo exp(2'0). (3) 

Varying ^ with respect to gives the equation of motion for the scalar field: 

= (^'^d^ijd'^^ - 2C,^e''^ + uj{i^)d^d^^^ + w(0)9'^V^^^^) • 

(4) 

We now specialize to the case of a spatially-flat Friedmann universe with 
coordinates and metric g^^ = diag(— 1, a^(t), a^(t), a^(i)). We wiU adopt 

a system of units in which G = 1 and c = 1. The modified Friedmann equation, 
obtained by varying ([2| with respect to the metric, is: 

^ - y (p™ (1 + |C|e-'^) + Pnb + Pre-^^ + p^, + Pa) (5) 



Here, an overdot indicates a derivative with respect to the comoving proper 
time t. In this equation, pA = ^ is a constant denoting the density of dark 
energy, p„f, oc is the contribution of non-baryonic matter that is assumed 
to form the cold dark matter, and p^ = t/)^. This equation is unchanged 
from the corresponding equation for BSBM, since any new contributions due 
to the V'-dependence of u will be of the form f^j^^, and we have = 0. 
Furthermore, ([5| can be approximated by a standard Friedmann equation at 
late times if the conditions 

|C|e-2'A < 1 and < p, (6) 

are satisfied then, where Pc represents the dominant form of matter in the 
universe during the epoch being considered. In this paper we will be mainly 
interested in the cold dark matter and A-dominated eras, in which pc = Pm 
and Pc — PA respectively. (Note that the first condition in ^ is automatically 
satisfied if -0 increases to infinity at late times. However, this is not a neces- 
sary requirement: for example, |^| could be small and could tend to a large 
constant.) 

We will call these two conditions the consistency conditions, and will use 
them to check that our late-time solutions are consistent with our assumption 
that the solution for the expansion scale factor, a{t), is unaffected by the small 
variations in tp that we are interested in modelling. In contrast, the evolution 
of ip is significantly affected by the cosmological expansion. 

The evolution equation Q for ^ becomes 

^ + 3i7V'+^^V? = -^e-^ (7) 
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and, after substituting uj'{tP)iJj — to, we can rewrite this as: 



,1/2^ y'^l 



where N = —2(^p„ia^ is a constant. In some sense, this equation is just the 
BSBM equation in disguise. This can be seen by defining a new time coordinate 
T with ^ = so that T = (uj~^''^dt. (Since w^i/^ > 0, this is a 

monotonic reparametrisation.) Then (18]) becomes 



dT \dT J 



(9) 



which is of the same form as the corresponding BSBM equation for -0, but with 
a different time parameter, T. However, this observation is of limited use for 
finding solutions, because we do not know the behaviour of ip{t) at the outset, 
and so we do not know the behaviour of uj{ijj{t)). 

Inspection of the conservation equation ([7| reveals an important feature of 
the evolution of ip (and hence a). In the case iV > 0, if ■0 vanishes anywhere, 
then ^ > at that turning point, so ip can be monotonic or have a minimum, but 
it cannot have a maximum. Similarly, in the case N < 0, ip can be monotonic 
or have a maximum, but it cannot have a minimum. Therefore, in either case, 
ip cannot oscillate, although some analyses [TS] of the linearised version of 
have been misled into deducing oscillatory behaviour for a. 

At present, uj{tp) is an arbitrary function. However, if the observed behaviour 
of a{z) and ipit) is monotonic, it should be possible (in principle) to obtain the 
form of uj{tlj) if our observational data are sufficiently accurate to allow us to fit a 
function for a(z), and hence determine V'(^)- This can be done as follows. First, 
note that equation Q can be re-written as a first-order differential equation for 
w(t): 

which has the general solution 

IN r 



Given an expression for ip{t), we can calculate the right-hand side of (10) in 
terms of t, and then write t = t^ijj) to obtain an expression for w (?/')• 

For example, in the matter era with a oc t^^^ and pm the general 

solution ( |Iol ) allows us to solve explicitly for uj{t), given a specified behaviour 

for ^{ty. 

uj{t) = 4^ / e-^^i:t^dt. 
(Vi2)2 J 

If, for example, ip = A\nt, where A> and A ^ 1, then 



-(^) = Kr3A) + t^-A(r3A)+^^ 
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On the other hand, if = \nt, then 

wft) = \{2N\nt + k)= e-^'''(2Nij + fc), 

where k is an arbitrary constant. It is easy to check that the consistency condi- 
tions (|6f hold in both cases (which imphes that the approximate cosmological 
solution a ^ t^/^ is valid) . 

III. A CHOICE OF ujitP) 

The coupling function a;('0) is a free function and its introduction generalizes 
the BSBM model, in which uj was a constant. This generalization introduces 
another way in which spatial inhomogeneity can be created in ip and a. We 
are interested here in situations where a;(V') remains non-negative and finite for 
all real values of i). We make the choice ^(V') oc e'^'^, where /i is a constant to 
be determined or constrained by observation. Other forms of w('0) either take 
negative values (e.g. ^('0) — "0): blow up at the origin (e.g. a;('0) = or are 
more complicated than the simple form above. Thus, in the remainder of this 
paper, we will restrict ourselves to considering the form 

L.(V) = c.oe'^^, (11) 

where wq and ^ are constants. Note that the BSBM theory is recovered in the 
/i — > limit. 

IV. SOLUTIONS IN THE MATTER ERA 



We need to investigate the cosmological evolution of tp in the epoch that 
includes the astronomical observations of quasars and low-redshift information 
on the possible evolution of a, so we begin our study in the matter-dominated 
era. We will assume that the universe contains a perfect 'dust' fluid, with zero 
pressure, and that the background scale factor is of the form a (t) = 00^2/3. We 
will need to check that the consistency conditions Q hold in order to ensure that 
we can ignore any corrections to the scale factor evolution from the evolution 
of ip{t). We will present exact solutions for certain special cases, but for the 
general behaviour we will use numerical investigations to provide insight into 
the evolutionary behaviour of a. 

From ([7|, the general evolution equation for ^{t) in a Friedmann background 
with expansion scale factor a{t) is 

V'+^V'^ + a-^^ 12 

2 a 

where w = loqc^^ and C — N/{LLiQa^). In the matter-dominated era, the scale 
factor evolves as a{t) — aot^^^, so this equation reduces to 

^+^^" + ^^=^2 ■ (13) 
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The consistency conditions are 

< land e'^'/'ipH^ < 1 as t cx). 



The equation ( 13 ) can be transformed into an Eniden-Fowler equation using the 
substitution = ^ In {Ot~^): 

= ti^ . 

2 

Note that 9{t) must be positive, and that always has the same sign. The 
consistency conditions become {9/t)~'^^^ <^ 1 and (j) — 9t^^^ ^ 1. The second 

condition imphes that 9{t) must grow at most as tint as t — > oo (and that if 
9{t) ^ tint then the constant of proportionahty should be very small). In 
applying these conditions, however, we note that in the standard cosmology the 
matter era lasts for a limited period of time before giving way to an accelerated 
expansion with p ~ —p. Thus, for example, if /i < 0, then a solution with 9 
growing asymptotically faster than t would be acceptable, as long as 9{t) ^ nt 
throughout the matter epoch (where k is a small constant). 
For two special cases, exact solutions can be found: 

Case SI: In the case /x = —2, the governing equation is ^ = This 
equation is equidimensional and has solutions of the form 9{t) — t^ ^ where k is 
a constant. If C = \, then 

9 = l3iVt + P2Vtl-at 

where we require that /32 > 0, or that (32 = and /3i > 0. This solution satisfies 
both consistency conditions. The corresponding solution for tp is 

V' = ^lnt-ln(/3i +/32 Int). 

In this (very) special case, ip therefore increases logarithmically at late times 
and the fine structure 'constant' evolves as 



exp(2-0) 



t t 



(/3i +/32lni)' (lnt)2' 
On the other hand, if C < j, then 

1- V1-4C l + v'l-4C 

9 = fiit 5 + p^t 5 

at late times, where, again, we require that ^2 > 0, or that /32 = and /3i > 0. 
In the former case, we require C > for the first consistency condition to hold; 
in the latter case, there is no condition on C. (As noted above, it is possible to 
have a consistent solution in the former case when C < 0, as long as P2 is very 
small.) The corresponding solution for -0 is 

^/^ = i (1 - Vr^ic) Int - In (/32 + Pit-^^^^ . 
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At late times, this solution tends to -0 = ^ (l — ^/]~^AC) Ini — (3it ^'-^ and 
a evolves as 

-fl-\/l-4C 

^ - - tl'VT^_ 



{a2 + ait-^/'^^y 

Note that, in this case, i(j decreases at late times if and only if C < 0. 
Finally, if C > j , we obtain 



6 ^t^ ( ;3iCos I- — hit j +/32sin I- — Ini 

In this case, physical solutions do not exist because 9{t) takes negative values. 

Case S2: In the case ^ = —4, the governing equation becomes 9 = —2Ct~^ 
with the exact solution 9 = —2Ctlnt + Pit + /32, where /3i and /32 are arbitrary 
constants. Since we require 9{t) to be positive (at least at late times), this 
means that a solution exists only if C < 0. The corresponding solution for ijj is 

^ = -iln ^-2Clni + + ^ 

and a ^ 1/ In t at late times. This solution decreases at late times, and therefore 
it is necessary to impose the requirement that — 2Clnt,„A ^ 1, where tmh is 
the time at which the matter era ends. 

We now consider general values of /i. We determine asymptotic forms for the 
solutions and illustrate the behaviour with numerical plots of ip against In(lnt)) 
for a wide range of initial conditions. Remarkably, there are just two simple 
asymptotes: 

Case 1: C > 0: Here, there are two cases depending on the value of /it. 
If /X < —2, there is a finite-time singularity at which ■0, and thus a, becomes 
unboundedly large. However, if /i > —2, all solutions tend to a common asymp- 
tote, similar to the behaviour found during the matter era in the standard BSBM 
scenario described in [8 (Fig. [T]). This result does not seem to depend on the 
initial conditions. The rate at which -0 increases along this asymptote can be 
made as gentle as desired by taking /i sufficiently large (Fig. [2|. 

We will now obtain an explicit formula for the asymptote by showing that, 
at late times, 

V-^ — *— ln((2 + ^)Clnt) 



and so 



(lni)2/(2+M) 

Note that this is consistent with the earlier BSBM result 0, ~^ ^ In (2Clnt) , 
which is recovered when /x — 0. 

In order to extract the asymptotic solution in finer detail we will proceed as 
in [8]. Using the substitution x = \n{t/ti), where ti is an arbitrary constant. 



we can transform ( 13 ) to the following equation: 



iP" + 7/;' + = Ce-(2+p)'^ (14) 
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where a prime denotes d/dx. We now write 



-0 = A\n{Bx) + ^ a„x"". 



Substituting this series into ( 14 ) gives 

oo 



n=2 

2 



+ ^ [Y.^ na„x-^-'j = |^ _ (2 + ^) a^x'^j . 

In order for the leading terms (in x^^) to match, we need —1 = —(2 + and 
thus A= ^.We also need A = CB-^^+aO^ which forces B = (2 + ^)C. This 
leads to 



1 , ^^^A2 



{2 + ^i)x V2 



A^ — A — aij x ^ + ((n — fiA)(n — l)a„_i — nan) x 



n=2 

2 



The successive terms of (a„)^]^ can now be chosen so as to make the left-hand 
side vanish. This requires 

4 + /i 1 (4 + m)^ , 

which leads to 

1 1 / . . _ \ 1 



exp ( -(2 + m) anx " ) 



(2 + ^)x (2 + ^)a; ^\ ' ±1 / (2 + ^^)2: 

as a; — > 00. Therefore, we have 



2 + Ai 

and 



CXD 

ln((2 + fi)Cx) + J2 



a = e^''' = [(2 + pL)C\n{t/ti)]'^ exp 



4 + /i 



(2 + /i)2 1n(i/ii) 



(15) 



to leading order. 

In order for this solution to hold during the matter era, the value of ti should 
be less than the value of t at the beginning of the matter-dominated era. 

Case 2: C < 0. The analysis again breaks down into two cases. If /i > —2, 
then at some finite future time, -0 becomes unboundedly negative, driving a to 
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zero; we ignore this case with a finite-time singularity. However, if /i < —2, 
the behaviour mimics the C > 0, /i > —2 behaviour, and solutions tend to a 
common asymptote 

^->-|-ln((2 + ^)Cln(<Ai)) 
2 + /i 

where, again, ti is an arbitrary constant. The analysis is similar to the above, 
but now -0 and a both decrease with time (Figs. |3]and|4]). 
Summary 

To summarize, there are three interesting conclusions to be drawn, assuming 
this form of ujiip) = luq exp(/i'0): 

(a) if a decreases during the late matter era, we must have fi < —2 and 
C < 0. Note that care must be taken in this case to ensure that the consistency 
conditions for the Friedmann equation are not violated, since e"^'^ = [(2 + 
/i)Clni]~'+f' is increasing and, if —4 < ^ < —2, then e^'^tp'^t'^ is increasing. 

(b) if a increases during the late matter era, we must have fi > —2 and 
C > 0. In this case both consistency conditions are automatically satisfied. 
This includes the original case of the BSBM theory, in which = 0. 

(c) it is possible to have special behaviour of a in the case — —2, and 
various conditions on C are required depending on whether a is decreasing or 
increasing. 

V. SOLUTIONS WHEN a{t) = aaf" 

We now consider the evolution of a{t) in situations where the cosmological 
scale factor exhibits general power-law behaviour a — aot"', where n > 0. The 
equation of motion for the scalar field becomes 

^ + + t,n (16) 

where, as before, uj = uj^e^^^ and C = iV/(a;oao). The substitution x = ln(t) 
then yields the equation 

i^" + (3n - 1)^-' + ^V-" = Ce-(2+A')^+(2-3n)x (^7) 

for tpix). Two separate cases then arise, depending on whether /i = — 2 or 
/x^-2. 

Case 1: In the case /i — —2, the evolution equation becomes 

V^-^^ + T^ = ^- (18) 

Let y — ip. Then 

C 3n 2 

y = ^y+y ■ 
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This is a Riccati equation, which can be solved by making the substitution 
y = —-. Then u satisfies 

.. 3n. C 

The solution of this differential equation is 

^ 35r^ I 2 - 3n / ^ 3^?^ I 2 -in ] 



first kind. There are two types of late-time behaviour. If n < | then 



where pi and /32 are arbitrary constants, and Jm{x) is a Bessel function of the 

2 

3 

/I -3n\ TT ^ 



2 cos 



2 




so, at late times, ?A — >■ const. + Ini and a — > i"^". As noted earlier, neither tp 
nor a can exhibit oscillatory behaviour because, when il) is zero in (18 1, ip can 
only have one sign (determined by the sign of C). Thus, ip cannot have both 
maxima and minima. 
If n > I , then 



2-3n / \ 2-3n / 

= Ai + ^2^^"^" 

where Ai and ^2 are new arbitrary constants. Therefore, at late times, 

V' -> const. - In (a + /3t^^^") ^ const. 

and Q!(i) — > constant. Hence, -0 increases as Int if n < |, but tends to a con- 
stant if n > |. The latter case includes the important scenario of a curvature- 
dominated Friedmann universe with n — \ and is indicative of the asymptot- 
ically constant behaviour of ip and a also to be expected in the dark-energy- 
dominated situation with n — ^ cx). We will discuss this scenario in Section 
VI. 

Case 2: In the case ^ ^ —2, we can define new variables by w = 
and z = i^' {x). This leads to the two-dimensional dynamical system 

w' = (2 — 3n)w — (2 + ii)zw 
z' = -^^z^ + {\-in)z + Cw. 
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which has three fixed points: {w,z) = (0,0), (0, ^(l-3n)) and ( ^^"^^"20(2+^)^2^""^^ ^ ^ 
We will call these points a, b and c respectively. 

Note that the fixed points a and b are asymptotic fixed points, since, by 
definition, w never actually reaches zero. On the other hand, point c exists only 
if (^~3")(3"M+i2w-4) ^ Q rpj^jg fixed point corresponds to the exact solution 



ih = In 

^ 2 + /i 



2-3n, 
M + 2 



(2-3n)(3n(^ + 4)-4)_ 

and this solution reduces to the one found in for the BSBM (/i — 0) case. We 
now investigate the stability of these points in order to understand the late-time 
behaviour. 

Critical point a. For (0, 0), the Jacobian matrix has the eigenvalues (2— 3n) 
and (1 — 3n). Therefore, if n 7^ | and n 7^ |, this critical point is hyperbolic. 
It is a saddle point if | < n < |, an unstable node if n < |, and a stable node 
if n > |. The case n < ^ is not of physical interest because it derives from 
background expansion with p > p. 

Critical point b. For ^0, ^(1 — 3n)^, the two eigenvalues of the Jacobian 
are (3n - 1) and |^ 3"A'+i2«-4 ^ ^ ^^^^^ Jacobian is T = 6n - 1 + 

^(3n — 1) and the determinant is 13 = (3n — 1) (sn + ^(3n — 1)^ . It is helpful 

to note that - 4D = {^j^i^n - 1) + 1^ ^ > if a* 7^ 4(1 - 3n). This means 

that, in general, the fixed point is either a saddle point (if D < 0) or a node (if 
T2 > 4D > 0). 

If n 7^ I and n ^ 3(1^+4) ' ^^^^ critical point is hyperbolic. Various types of 
behaviour can occur depending on the values of n and fi: 

• 71 > I . Then this point is a saddle if 3n < ^ and an unstable node if 



3n > 



4-12n 



• n < Then this point is a saddle if 3n > ^—^^ and a stable node if 
3„ < i^. 

Critical point c. In this case, the Jacobian matrix is 

(3n-2)(3Kp+12K-4) \ 
2(2+p)C I 

2+M / 

SO T = and D = (2 3n)(3ri.|j+i2n-4) ^ bcforc, there are several cases. 

Note that we do not distinguish here between whether the point is a node or 
a focus, since our primary interest is whether the critical point is stable or 
unstable. 

• If n > I and /i < 1^ — 4, point c is a saddle point. 

If 71 > I and ^ — 4 < /i < —2, point c is an unstable node/focus. 
If 77 > I and /i > —2, point c is a saddle point. 
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• If n < I and /i < —2, point c is a stable node/focus. 

If 71 < I and —2 < jj < ~ 4, point c is a saddle point. 

If n < I and ii > ^ — 4, then point c is a stable node/focus if ^ > 2 — 6n, 
and is an unstable node/focus ii fi < 2 — 6n. (Note that this means that 
if n < ^ and /i > — 4, point c must be a stable node/focus.) 

We can use these results to obtain a general result for the case n > | 
by noting that the only stable critical point is at {w,z) — (0,0). Therefore, 
assuming that z = tip does not diverge to infinity, we are likely to obtain the 
behaviour tp — >■ const, at late times in accelerating universes. 

VI. SOLUTIONS IN AN ACCELERATING UNIVERSE 

The data from the Oklo natural reactor {z « 0.14) |20j, the analysis of 
meteorites {z « 0.45) [H] and laboratory experiments [H] provide very tight 
bounds on the allowed variation of alpha at low redshifts. We will show that our 
generalization of the BSBM theory is still consistent with negligible variation in 
a during the accelerating era of the universe, although we do not introduce the 
added complexity of relating local terrestrial observations of varying constants 
to their global cosmological evolution (see Barrow and Shaw [53], for 
a detailed discussion of this problem). 

The cases of power-law expansion with n > 2/3 discussed in the last section 
cover the situation of an accelerating universe dominated by a fluid with equa- 
tion of state —p < p < — /o/3. The remaining situation of special interest is the 
particular case of a A-dominated era, with p = — p, in which the scale factor 
evolves according to the dc Sitter form 



where M = •v/A/3. The evolution equation (12) for ip becomes 



^ + ^^^ + 3M^=^— ^ (19) 
where C = N/{ujoaQ), and the consistency conditions ^ are e~'^^ <^ 1 and 

Numerical investigation shows that there are two types of solution. One 
type tends to a constant at late times, in general increasing before it does so. 
The value of ■0 along the eventual constant asymptote depends on the initial 
conditions. A notable feature of Fig. [5] is that the behaviour of during the 
initial stages of the evolution can be complicated, but all solutions seem to 
become constant for t > t* , where t* is some fixed 'cooling-off' time. Therefore, 
if the matter era lasts past time t* , we would expect negligible variation in ip 
once the cosmological constant begins to dominate. 



We might expect this behaviour from ( 19 ) if -0 is not too large, because then 
we have, approximately, 

^ = -3Mip. 
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Since 3M > 0, i/' is driven exponentially quickly to zero, and ip tends to a 
constant. 

The second case occurs when ip is large, so that the ^tp'^ term comes into 
play. We then have 

= - (IV- + 3m) V', (20) 

and there are two equilibria for ip: at -0 = and ip = —QM/ ^. Regardless of the 
sign of fi, the fixed point a.t ip = is stable and the fixed point at tp = — ^ is 
unstable. Therefore, depending on the initial conditions, there are two generic 
behaviours for tJj: it may either tend to ±oo exponentially quickly, or tend to 
zero. Assuming that finite-time blow-up does not occur, ^p should tend to zero 
at late times, so ^ tends to a constant as found above. Otherwise, we obtain a 
second class of diverging solutions. 

We can also perform an analysis similar to the one performed in Section V 
for the general power-law case, but this time we will work directly with equation 



(19). We begin by defining the variables w ^ e 3A/t ^^j^j ^ ^ rpj^jg 



leads to the autonomous dynamical system 

w — —3Mw — (2 i.i)zw, 
i = -^z'^ -3Mz + Cw. 

Again, this system has three fixed points, at (Wjz) — (0,0),(0, — ^) and 

(~ ^2C{2+Ij.)^ ' ~ j^A' ) • before, we will name these points a, b and c respec- 
tively. Again, points a and b are asymptotic fixed points, and point c only exists 
if (4 -I- /i)/C < 0. The analysis of the stability of these fixed points now follows: 
Critical point a. At (0,0), the Jacobian has the two eigenvalues — 3Af. 
Therefore, the origin is always a stable node. 

Critical point b. At ^0, ^ the Jacobian has eigenvalues 3M ^1 + 

and 3Af. The trace and determinant are therefore given by T = 6M (^^^ + 1 

and D = 9M^ + 1^ respectively, and we have - 4L) = 144 > 0. Hence, if 
—4 < fi < 0, this critical point is a saddle. Otherwise, ii fi > ot fi < —4, it is 
an unstable node. 

Critical point c. If this fixed point exists, the corresponding Jacobian 
matrix has determinant D = —^^w?4p-v^ and trace T ~ — tt^- Therefore, if 

212+fj,) 2+11 ' 

—4 < n < —2, this point is an unstable node, and if cither ji > —2 or < —4, 
it is a saddle point. 

Note that, for any value of M, the origin (0, 0) is the only stable node. 
Therefore, in an accelerating era, any bounded late-time solution will satisfy 
Ip — >■ const., and we have a form of the cosmic no-hair theorem for the evolution 
of ■0 and a in the de Sitter background metric. 

VII. SUMMARY AND THE RELATION TO OBSERVATIONS 
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We are now in a position to summarize the results of the previous sections, 
having determined the late-time behaviors for both power-law and exponential 
evolution of the scale factor a{t). For convenience, the different outcomes are 
displayed in Table 4.1. 



Value of 


Scale factor evolution a{t) oc t" 


a{t) oc e^^' 


5<"<l 




n > 1 


H>0 


2(2-3ii) 

Of -)• t 2 + ^ if C > 


Q [(2 + ^)Cln{t/ti)]^ ifOO 


const 


const 


IJ. = (BSBM) 


a t^-^" if C > 


a 2Chi{t/ti) if C > 


const 


const 


-2 < ft<0 


2(2-371) 
Of t 2 + ^ if C > 


a [(2 + ^)Cln(t/ti)]^ if C> 


const 


const 


^< = -2 


a^t^" if C > 


various cases, depending on C 


const 


const 


fj.<-2 


2(2-371) 

a i if C < 


[(2 + ^)Cln(t/ti)]^ if C<0 


const 


const 



Table 1: This table lists the late-time behaviours of a{t) for various choices 
of fj, and possible evolutions of the scale factor a(t), of both power-law and 
exponential forms. In the power- law case with n > 2/3, and in the exponential 
case, the fine structure 'constant' a{t) tends to a constant value at late times. 
In other cases, a{t) tracks an increasing asymptote if /i > —2, and a decreasing 
asymptote if fi < —2. 

We can now attempt to constrain the parameters of the theory by using 
observational data [26j . Since most of the observational evidence arises from 
the late matter-dominated era, we will use the asymptotic behaviour obtained 
in earlier sections as the theoretical basis for comparison. 

In the matter-dominated era, the redshift-time relation 

ajh) ^ 1 + Z2 ^ /iiy^' 
holds. From our results in Section IV above, solutions tend to the asymptote 

where ti is an arbitrary constant (corresponding to redshift zi, say). Then, 
using Q, we can express the relative shift in the value of the fine-structure 
constant expected between a redshift z > z\ ^ 0.5 and the present by 

Aa a{z) — ao a{z) — aA 

/ ln(l + zi) - ln(l + ^) \ ^ 
[\nil + z,)-\nil + ZA)) 
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where aA = Q^o (to an excellent approximation) is the value of a at the redshift 
signifying the onset of A-domination and an accelerating universe. 

Thus, defining the new parameters i/ = 2/(2 + /i) and Z = ln(l + zi), we 
obtain the following approximate form for the evolution of Aa/a during the 
matter era: 

Aa ( Z - ln(l + z) 



- 1. 



(21) 



,Z-ln(l + ZA), 

The only constraint on Z is that it should correspond to a redshift before dust 
domination, i.e., Z > 10. Taking za = 0.5, sampling values of Z between 
0.1 and 1000, and values of i' between —100 and 100, indicates that the best-fit 
parameters generally satisfy the approximate relation v ^ 2.5Z, or {2 + fj,) ln(l-|- 
Zi) ^ 0.8. The associated value is approximately 168.75, with 143 degrees 
of freedom; this indicates a fit that is reasonable but not particularly strong. 



Formula (21) generalises the analytic approximation for Aa/a obtained in the 
BSBM theory (^i = 0). 

In the above, we have assumed a sharp transition between dust domination 
and de Sitter expansion at z = za, and assumed that ijj and a are constant 
when z < z\. A more accurate approximation that evolves smoothly between 
the dust and A-dominated eras could be obtained by assuming a background 



cosmological evolution in (12 1 where 



3 SnGpoa^ . 2 
a = : smh 



A 



/3A 



(22) 



This form of a{t) solves the Friedmann equation 

d2 SttG a 



for a universe containing dust and a cosmological constant. In this model. 



the fractional densities are Qm = sech^ (^^^^) ^^"^ ~ tanh^ 
spectively. Using the approximate values 17 ao = 0-73 and 
and Ho = 74.2kms~^Mpc"^ [Ml, we obtain A 



/3A. 



t] re- 



0.27 



1.266340 X 10-35s-2 and 



SirGpo = O.Sli^o^ = 4.684 x lO^-^^s-^ We have data for 



Aa 

a 



-24,0 



/-I 



and we are interested in the evolution of f{z) = e 



2-41 



The evolution equation for ip{t), (13), can be transformed into an evolution 
equation for f{z): 



[SnGpoil + zf +A] (1 
A8TrGDpoil + z)^f-^, 



4/ 



/ 



- [47rGpo(l 



2A] /' 
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where D — —-r-r-^ — ■ This allows us to determine the behaviour of f for various 
initial conditions /(O) and /'(O). Note that /(O) = e^'^f > 0, and, since the 
present-day variation of a is neghgible, we have assumed that /'(O) — 0. Thus, 
our system has three arbitrary parameters: the value of /(O), the value of 
/i, and the value of D. The contour plots in Fig. [6] illustrate the regions of 
parameter space for which there is a good fit to the observed data. The best 
chi-squared statistic obtained using this model is approximately 169 (attained 
when /(O) < 1), which still indicates a fit that is not very strong. 

It is important to note that these investigations of the time variation of a 
may not be conclusive because recent observations of a are consistent with a 
spatial variation, and the latter will need to be accounted for before an accurate 
estimate of jj, (or an accurate test of whether /i ^ 0) can be obtained. If we 
are interested in modelling possible spatial variations in a at the redshifts of 
quasar absorption spectra, the BSBM theory would require inhomogeneity in 
C,/uj. Since uj is constant in that theory, it would be necessary for C, to vary in 
space. This means that the identity of the dominant form of cold dark matter 
would vary throughout space, or has a strange character that requires spatial 
variation in its electric and magnetic composition. This is not appealing in 
the absence of other evidence because the time evolutions of the densities of 
these different varieties of pressureless matter will be the same even if they have 
different C, values, and, for one value of C to be replaced by another, some type 
of decay would need to occur. However, in the generalized theory explored here, 
inhomogeneity in a can be caused by spatial variation in (^/u!{ip) via ^^{x), and 
this is not unnatural. This situation will be explored in greater detail elsewhere 
by extending the results of 29J and |30j . 

VIII. CONCLUSIONS 

We have shown that the BSBM theory of varying alpha naturally generalises 
to cases in which the scalar field coupling a;(V') > has a functional dependence 
on the scalar field ip. We have explored the features of the specific case where this 
dependence is exponential, w('0) c>c e^"^ . The BSBM solution, which is recovered 
when /i = 0, may now be viewed as just one of a range of possible solutions with 
jjL > —2 which all exhibit similar asymptotic cosmological behaviour in dust or 
dark-energy-dominated eras at late times. In all these cases, a increases slowly 
with time, or (in the dark-energy-dominated cases) increases and then tends to 
a constant-valued asymptote. In the case of a de Sitter background, this gives 
an extension of the cosmic no-hair behaviour familiar from studies of general 
relativistic cosmology with constant a. 

We have also found different possible behaviors in the cases < —2. In these 
cases, it is possible for ip (and a) to decrease with time, and the consistency 
conditions must be checked numerically in order to ensure that our approximate 
solution to the Friedmann equation is not changed by the growing gravitational 
effects of decreasing tp. 

Finally, it is important to note the limitations of this analysis. In this paper, 
we have confined attention to theories where the late-time behaviour of a{t) is 
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unaffected by the varying scalar field. However, if this was not the case, then 
new possibilities might arise. We have also highlighted the need for an analysis 
of the effects of spatial variation of ip and a in the post-recombination era; this 
work will be presented elsewhere. 

ACKNOWLEDGMENTS 

S. Lip would like to thank the Gates Cambridge Trust for its support. We 
also thank Anne Davis, Joao Magueijo, Michael Murphy and John Webb for 
discussions and Sergei Levshakov for helpful comments.. 



References 

[1] J. K. Webb, V. V. Flambaum, C.W. Churchill, M. J. Drinkwater and J.D. 
Barrow, Phys. Rev. Lett. 82, 884 (1999). 

[2] M. T. Murphy, J. K. Webb, V. V. Flambaum, V. A. Dzuba, C.W. Churchill, 
J. X. Prochaska, J.D. Barrow and A.M. Wolfe, Mon. Not. Roy. Astron. Soc. 
327, 1208 (2001). 

[3] J. K.Webb, M. Murphy, V. Flambaum, V. Dzuba, J.D. Barrow, C. 
Churchill, J. Prochaska, and A. Wolfe, Phys. Rev. Lett. 87, 091301 (2001). 

[4] M.T. Murphy, V.V. Flambaum, J.K. Webb, V.V. Dzuba, J.X. Prochaska 
and A.M. Wolfe, Lect. Notes Phys. 648, 131 (2004). 

[5] J.K. Webb, J.A. King, M.T. Murphy, V.V. Flambaum, R.F. Carswell and 
M.B. Bainbridge, arXiv:1008.3907. 

[6] J.D. Bekenstein, Phys. Rev. D 25, 1527 (1982). 

[7] H.B. Sandvik, J.D. Barrow and J. Magueijo, Phys. Rev. Lett. 88, 031302 
(2002). 

[8] J.D. Barrow, H.B. Sandvik and J. Magueijo, Phys. Rev. D 65, 063504 
(2002). 

[9] C. Brans and R.H. Dicke, Phys. Rev. 124, 924 (1961). 
[10] J.C. Berengut et al, arXiv:astro-ph/1011.4136. 

[11] S.A. Levshakov, P. Molaro, S. Lopez, S. D'Odorico, M. Centurion, P. Boni- 
facio, I.I. Agafonova and D. Reimers, Astron. Astrophys. 466, 1077 (2007). 

[12] P. Molaro, D. Reimers, I.I. Agafonova and S.A. Levshakov, Eur. Phys. J. 

Special Topics 163, 173 (2008). 

[13] C. Martins, arXiv:astro-ph/1011.6591. 

[14] J.-P. Uzan, Living Reviews Relativity 14, 2 (2011). 



18 



[15] R. Srianand, H. Chand, P. Petitjean and B. Aracil, Phys. Rev. Lett. 92, 
121302 (2004). 

[16] M.T. Murphy, J.K. Webb and V.V. Flambaum, Phys. Rev. Lett. 99, 239001 
(2007). 

[17] J.D. Barrow, H.B. Sandvik and J. Magueijo, Phys.Rev. D 65, 123501 
(2002). 

[18] E. J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 
(2006). 

[19] J.D. Barrow and D.F. Mota, Class. Quant. Grav. 19, 6197 (2002). 

[20] T. Damour and F. Dyson, Nucl. Phys. B 480, 37 (1996). 

[21] K.A. Olive, M. Pospelov, Y.-Z. Qian, A. Coc, M. Casse and E. Vangioni- 
Flam, Phys. Rev. D 66, 045022 (2002). 

[22] S. Blatt et al, Phys. Rev. Lett. 100, 140801 (2008); T.M. Forticr et al, 
Phys. Rev. Lett. 98, 070801 (2007); E. Peik et al, Phys. Rev. Lett. 93, 
170801 (2004); E. Peik et al, Proceedings of the 11th Marcel Grossmann 
Meeting, Berlin, 2006, World Scientific, Singapore (2008); M. Fischer et 
al, Phys. Rev. Lett. 92, 230802 (2004); A. Cingoz et al, Phys. Rev. Lett. 
98, 040801 (2007); T. Rosenband et al. Science 319, 1808 (2008). 

[23] D.J. Shaw and J.D. Barrow, Phys. Rev. D 73, 123505 (2006). 

[24] D.J. Shaw and J.D. Barrow, Phys. Rev. D 73, 123506 (2006). 

[25] J.D. Barrow and D.J. Shaw, Phys. Rev. D 78, 067304 (2008). 

[26] M.T. Murphy, J.K. Webb and V.V. Flambaum, arXiv:0911.4512. 

[27] E. Komatsu et al (WMAP collaboration), Astrophys. J. Suppl. 192, 18 

(2011). 

[28] A.G. Riess et al, Astrophys. J. 699, 539 (2009). 

[29] J.D. Barrow, J. Magueijo and H.B. Sandvik, Phys. Rev. D 66, 043515 
(2002). 

[30] J.D. Barrow and D.F. Mota, Class. Quant. Grav. 20, 2045 (2003). 

[31] J. Magueijo, J.D. Barrow and H.B. Sandvik, Phys. Lett. B 549, 284 (2002). 



19 





Figure 1: A plot of four solutions of equation (13 1 for the various choices 
'ijj{lnx — 0) £ { — 1,-0.5,0,0.5} and d^/j / d{\n x)\in x=o = 0. Here, x is defined 
as ln(i/ti), where ti is an arbitrary constant, and we have chosen (arbitrarily) 
the parameters n — 2 and C — 0.1. Notice that, in this case (/i > —2), all the 
solutions eventually converge to a common increasing asymptote, independently 
of the initial conditions. 
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Figure 2: As in Figure [T| but with /i = 6. All solutions still converge to a 
common increasing asymptote at late times. However, the slope of the asymp- 
tote is more gentle than that in Figure [l] because the value of /i has increased. 
The slope of the asymptote can be made as gentle as desired by taking to be 
sufficiently large. 
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Figure 3: A plot of four solutions of equation ( 13 1 for the various choices 
ip(hix = 0) e {—1,-0.5,0,0.5} and dip/d{\nx)\\nx=o = 0, where, again, x 
is defined as ln(t/ti). The parameter values are now /i = —4 and C = —0.1. In 
this case (/i < —2), the solutions eventually converge to a common decreasing 
asymptote, and this behaviour is insensitive to the initial conditions. 
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Figure 4: As in FigurejSj but with /i — —6. All solutions still converge to a com- 
mon decreasing asymptote at late times. However, the slope of the asymptote is 
more gentle than that in Figure [S] (because /i is now more negative). Again, this 
slope can be made as gentle as desired by taking a sufficiently large negative /i. 
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Figure 5: As in Figure [T] but these solutions now occur within an accelerating 
universe. The parameter values are now fjL — 2, C = 0.1 and M = 0.01. 
Regardless of the initial conditions, the solutions all tend to constant values after 
a short initial time. This illustrates how the time variation of a is suppressed 
in a A-dominated universe. 
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Plot of mu vs □ with f[0) = 0.001 plot of mu vs D with f[0) = 0.01 




Figure 6: These contour plots display the value of the statistic (according 
to the legend shown in the right-hand color bars) for various choices of D, fj, 
and in a model that exhibits a smooth transition between the dark-matter 
and dark-energy dominated eras. Each picture corresponds to a different value 
of /(O) = e^'''^''^ ranging from 10^^ to 10^. The vertical and horizontal axes 
represent the parameters ^ and I?, respectively. Darker regions in the plots 
indicate parameter choices that give rise to a lower statistic, and therefore 
correspond to a better fit. 



25 



